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1. Introduction 

Exact solutions of quantum field equations of motion provide an effective tool in 
studies of different phenomena of particle interactions in high energy physics. They 
supply with particular applications in solving problems of charged particles motion 
in electromagnetic fields of terrestrial experimental devices, as well as in astrophysics 
and cosmology. Exact solutions were first applied in quantum electrodynamics for 
development the quantum theory of the synchrotron radiation, i.e. for studies of motion 
and radiation of the electron in a magnetic field (see, for instance, [I]), and also for 
studies of the electrodynamics and weak interaction in different other configurations of 
external electromagnetic fields [2]. This method is based on the Furry representation [3] 
in quantum electrodynamics (for more detailed discussion on this item see [3]) widely 
used for description of particles interactions in the presence of external electromagnetic 
fields. Recently it has been shown, that the method of exact solutions can be also 
applied for the problem of neutrinos and electron motion in presence of dense matter 
(see [1] for a review on this topic). Most pronouncedly this possibility was pointed out 
in [HI [6] where the exact solution for the modified Dirac equation for a neutrino moving 
in matter was derived and discussed in details. The corresponding exact solution for 
an electron moving in matter was obtained in [7], the problem of neutrino propagation 
in transversally moving matter was first solved in [8], and in our recent paper [9] we 
considered neutrino propagation in a rotating matter accounting for the effect of nonzero 
neutrino mass. 

Here we further develop the method of exact solutions for the problem of charged 
leptons propagating in matter and strong magnetic fields. The paper is organized as 
follows. In Section 2 we discuss the general form of the modified Dirac equation for 
an electron moving in matter and magnetic field, and the corresponding spin operator 
is derived in Section 3. In Sections 4-6 the exact wave functions and energy spectrum 
are obtained. The general problem of eigenvectors and eigenvalues of a given class of 
Hamiltonians are discussed. In conclusion (Section 8) we examine the exact solutions 
obtained in Sections 2-6 for the particular case of the charged particle motion in a 
rotating neutron star with account for matter and magnetic field effects. 

2. Modified Dirac Equation for electron moving in matter and magnetic 
field 

We consider an electron propagating in nonmoving magnetized medium composed of 
neutrons and suppose, that magnetic field is homogeneous and constant. This can be 
regarded as the first approach to modelling of an electron propagation inside a neutron 
star. For distinctness we consider here the case of an electron, whereas generalization 
for other charged particles is just straightforward. We start with the modified Dirac 
equation for the electron wave function exactly accounting for the electron interaction 
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with matter in the absence of magnetic field [5] (see also [5]): 

{l^ + ~7„(1 - 4 sin 2 6 W + 7 5 )r - m}*(z) = 0. (2.1) 
This is the most general form of the equation for the electron wave function in which the 
effective potential = \{1 — 4 sin 2 6w + 7 5 )/^ includes both the neutral and charged 
current interactions of the electron with the background particles, and which can also 
account for effects of matter motion and polarization. 

In order to include also an effect of an external electromagnetic field, we replace in 
eq. (12. ip the electron momentum p^ by the "extended" momentum: p M — > p^ + CqA 11 , 

{ 7 „(p" + e A») + 1 7m (1 - 4 sin 2 9 W + 7 5 )/^ - m}*(x) = 0, (2.2) 

where eo is a module of the electron charge. 

Note, that in general case it is not a trivial task to find solutions of this equation. In 
what follows we consider the particular case of constant magnetic field and nonmoving 
uniform matter, so that for the electromagnetic field and effective matter potential we 
obtain 

^ = (0,— ,np0), r = -Gn(l, 0,0,0), (2.3) 

where G = n is matter number density. We rewrite equation (12. 2 j) in the 

Hamiltonian form and get 

ij^tf(x) = (2.4a) 

# = 7°7(P + e A) + m 7 ° + -(1 - 4 sin 2 6W + 7 5 )Gn (2.46) 

where A = (— -^,^,0). Using the chiral representation of the 7 -matrices we obtain 
the Hamiltonian in block matrix form: 

H=( -<r'(P + e oA) + Gn{l-2sm 2 6 w ) m \ 

I m o-'(p + e A)-2Gnsm 2 6 w ) [ ' ' 

where cr' are Pauli matrices. This form of the Hamiltonian makes quite transparent the 
solution describing spin properties of the electron. 

3. Spin operator 

It is obvious from ( 12. 5p . that the longitudinal polarization operator T° = ^er(p + e^A) 
[1] , where cr = [ ^ \ are the Dirac matrices, can be written in the form 



fo = ^'(P + e A) 

i-cr'(p + e A) 



(3.1) 



and commutes with the Hamiltonian, [T°, H] = 0. Therefore for any of its eigenvectors 
the Hamiltonian can be presented as the matrix: 

„ . -mT° -2Gn sin 2 6 W + Gn m \ 

! m mT° - 2Gnsin 2 ^ j ' ( j 
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where T° is one of the eigenvalues of the spin operator T°. Note, that the matrix (13. 2p 
is still 4 ® 4 one, and each its element is actually a product of the number by 2 £g> 2 
unit matrix. In turn, the transverse polarization operator does not commute with the 
Hamiltonian (12. 5p because of the matter term Gn that breaks block symmetry of the 
Hamiltonian. 



4. Energy spectrum of electron in matter and magnetic field 

In order to find the electron energy spectrum po in the matter and constant magnetic 
field, H^/ = Po^/, we should solve the equation 



0, (4.1) 



-mT° + Gn — po m 

m mT° — p 

where po = Po + 2Gnsin 2 9w- The solutions can be written in the form 



Po = 5p - 2Gn sin 2 9 W + eJ (mT° - + m 2 , (4.2) 

where e = ±1 is the "sign" of the energy. 

It is significant to note an interesting feature of the electron energy spectrum in the 
magnetized matter following from (14. 2p . It is well known, that the energy spectrum of 
the electron in the magnetic field is degenerated in respect of spin quantum number (each 
electron Landau energy level in the magnetic field corresponds to both spin orientations). 
The presence of the matter (of any non- vanishing density n ^ 0) removes the degeneracy. 
This phenomenon can be attributed to the parity violation in weak interactions. 

Let us emphasize one important relation between po and To that immediately follows 
from the spectrum (14. 2 p : 

(po _ 9?L + 2Gnsin 2 9 w f = (mT° - ^) 2 + m 2 , (4.3) 

where T° is one of the eigenvalues of the spin operator T°. We use this relation in 
Section 7. 



5. The electron wave functions 



Note that we can considerably simplify the problem of finding wave functions if take 
into account some obvious facts. The solution of the equation (I2.4ap due to symmetries 
can be sought in the form 

-i Po t+- m z fo(x,v) 



V(t,x,y,z) = e 



(5.1) 
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Substituting ( 15.1 j) into (I2.4al) we arrive at a system of linear equations for the electron 
wave function components: 



d_ .d_ 

dx dy 



dx dy 
e B 



(On - /;;.)''! 4- i <J ( - i^- ) + ~^2~( x ~ iy) } 1 '-2 ••• i - ( r >.2«) 



x + iyH -01 + (P3 + Gn)^ 2 + m0> 4 = Po^2, (5.26) 



a . 8 



dx dy 



e B, 



imp! + p 3 ip 3 + i < ^ + — -(x - iy) } ip A = p ^ 3} 



d . d 



e B 



rmjj 2 + i | + i-^J —{x + iy) j ip 3 - £>30>4 = £»o04, 

In the polar coordinates x + iy = re 1 ^, x — iy = re~ lr ^ one obtains 



JL _i_ jJL — e i0 [ o_ _|_ j__c<_ > _o _ e -i0 i _o_ _ j__a_ | 

f)^r. d r u \ dr r d(h I ' fir. fin \ fir r firh I ' 



dx dy 



d , i d 



dr r d(f> 



d ■ d 



dx dy 



d i d 



dr r d<j> 



and the system of equations (15 .2 aft - (15.2dj) reads now 



-/);. : (-///)( io l! ' <j - X -^- + ^ i ' 2 - ///;•,, /">„< ' j . 
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d i d e B , 

+ -— —r > ip! + {p 3 + Gn)i) 2 + mip 4 = p ip 2 , 



dr r d(J) 2 
mipi + P3"03 + ie" 1 ^ 

d i d e B 



rmjj2 + ie 1 ^ 



dr ^ r 



d i d e B 

^ m. + ~^r r f ^ 4 =po^3, 

or r do 2 



r > ^3 -P3^4 = PO04- 



(5.2c) 
(5.2d) 

(5.3) 

(5.4a) 
(5.46) 
(5.4c) 
(5 Ad) 



It is possible to show that the operator of the total momentum J z = L z + S z , where 
L z = —i-§^i S z = |<73, commutes with the Hamiltonian. Therefore the solutions can be 
written in the form 



02 
0>3 



ix 2 (^)e 1 ^ 
\ ix4(r)e ll<t ' J 



(5.5) 



The solutions are the eigenvectors of the total momentum operator J z with the 
corresponding eigenvalues I — |. After substitution of f)5.5p . the system (15 .4 aft - fl5.4aT) 
turns into: 



(p 3 - Gn)xi - 
d l-l e B 



dr r 2 

mxi + P3X3 + 

d l-l e B 



d / e B 

+ ~ + — r J X2 + "2X3 = P0X1, 

r I Xi + (P3 + Cn)x 2 + mx 4 = P0X2, 



d Z e 5 



mx 2 



dr r 



r X3 -P3X4 =PoX4- 



(5.6a) 
(5.66) 
(5.6c) 
(5.6d) 
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Now we define the "increasing" and "decreasing" operators 
d l-l e B 
dr 



R- 



r 2 

and get the Hamiltonian in the following form, 

f-p 3 + Gn -R- m \ 

R + p 3 + Gn m 

m pz R~ 

\ m -R + -psj 



d I e B 

- — I 1 r. 

dr r 2 



(5.7) 



H 



2Gn sin 2 9wl, 



(5. 



where I is the unit matrix. Note that the derived form of the Hamiltonian is spectacular 
transparent that significantly simplify the problem of getting the explicit form of 
eigenvalues of the Hamiltonian ( 12. 4 all . For forthcoming applications we take into 
consideration properties of operators R + and R~: 

'e B 



R+£ l s -i 



R C\ 



e B 2 
r 



'2eoB{s + l)C l a 



e Q B 



2e B(s + l)C 



i-i 



e B 



(5.9a) 
(5.96) 



where C l s are the Laguerre functions p]. 

The solution of system ( 15.6 ah - ( 15.6 all (the eigenvector of the Hamiltonian ( 15.8 
can be written in the form 

I Xx \ ( C.C 1 - 1 ( e -fr 2 ) \ 



X'2 
X3 
\X4/ 



C 2 L\ { e -fr 2 ) 



l-l (e B r 2 



(5.10) 



(-2 

Now we can get the eigenvalues of the Hamiltonian (the energy spectrum) and of 
the spin operator T°, 



Gn 



2Gn sin 2 6 W + e 




p\ + 2e B(l + s) +m 2 , e±l, (5.11) 



■m 



p 2 3 + 2e B(l + s), s' = ±l. 



(5.12) 



It is easy to see, that the spectrum (14. 2 p obtained above is in agreement with expressions 
( 15.1 ip and ( I5.12p . From this energy spectrum, it is straightforward that the well-known 
energy spectrum in magnetic field (the Landau levels) is modified by interaction of the 
electron with matter. However the radius of the classical orbits corresponding to a 
certain level (15. lip doesn't depend on the matter density: 



(R) 



v]>+ r v|> dr 



2N 



(5.13) 



To conclude this section we would like to note, that the effect of electron trapping on 
circular orbits in magnetized matter exists, and this can be important for astrophysical 
applications. 
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6. General problem of eigenvectors and eigenvalues of a given class of 
Hamiltonians 



We argue in this Section, that the problem considered above is a particular case of 
more general problem, which concludes in a retrieval of eigenvectors and eigenvalues of 
Hamiltonians with a given general structure. The examples of these particular cases 
can be found in papers [U |8j [9]. Taking into account the increasing interest to such 
Hamiltonians, we expect that their profound consideration would be important for 
applications. 

Theorem . Let us take into account the following conditions: 

1) % is a Hilbert space with the basis of eigenvectors ip{n} of given quantum problem 
for a Dirac equation ( a matter and magnetic field are included); 

2) in this space, the increasing and decreasing operators are determined 

aip{ n } = /-(W)VV-i}' ( 6 - la ) 
a + ip{ n -i} = f+{{n})ip{n}, (6.16) 

where /_({n}) and f + ({n}) are known functions that depend on the set of quantum 

numbers {n}; 

3) the Hamiltonian has one of the following explicit structures 
/ m v a \ ( - 



H 



m 





P 


a 





m 


a + 


~P 


P 


a 


—m 





a+ 


~P 





—m 



or 



H 



-p 


a 


m 


o\ 


a+ 


P 





m 


in 





P 


—a 





m 


-a + 


-p) 



(6.2) 



where m is the mass of a particle, p is (often a third) component of momentum, although 
here we could use other symbols for demonstration the general mathematical structure. 
Note that the theorem is also valid in other case when a certain class of replacements 
within the matrixes (16. 2 p are made, for example, a -H- a + , a — > ia, a + — > — ia + , p -H- — p, 
m -H- — m, etc. Nevertheless, the disposition of operators a, a + , and zeros in the matrixes 
important here. Note that there should be only one operator (a or a + ) and only one 
zero in any string and any column. Remarkably, for any string or column the operators 
and zeros occupy only even or only odd positions. 

In that case the following statements can be proved: 

1) the eigenvector of such a Hamiltonian has the following form 

( CiVv-i}\ 

(6.3) 



C 2 ?/>{„} 

\ C 4 ^{n} J 

and the equation for the spectrum is (for the left structure in f)6.2p ) 

m — E p /_({n}) 

m — E / + ({n}) — p 
p f-({ n }) —m — E 

/ + ({n}) —p —m — E 







(6.4) 
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The equation for the right structure in (16. 2p can be written in analogous way. The 
spectrum can be obtained explicitly in the form 

E = ey/rn?+p* + f + ({n})f_.({n}), e = ±1; (6.5) 

2) the spin operator (one of the possible variants) can be constructed from the 
Hamiltonian. It is a matrix of operators, zeros, integrals of motion and other parameters 
(m). A structure of the spin operator corresponding to the Hamiltonians in ( 16. 2p can 
take one of the following variants: 

/* 



S 



/* 




* 

1+ 





* 
a 



a+ 
* 




* 


* 



or 



S 



* 




a 
* 

* 



* 




(6.6) 



where all free positions * should be filled in using the main condition [S, H] = 0. Of 
course, some refinement of structure of spin operator (or any of its blocks) is possible, 
for example, a — > ca, a + — > c*a + , where c is a complex number, |c| = 1. Note that the 
main part of known spin operators (see, for example [TJ) has such a structure. 
The proof of the theorem 

1. Let us substitude the solution (16. 3p into equation Hip = Eip where H is given by 
one of the explicit structures (16.21) . The homogenious system of linear equations is now 
obtained. The equation (16. 4ft is a criterion of its nonzero solution, and it is the equation 
for energy spectrum. To obtain the spectrum (16. 5p we should solve this equation, what 
is trivial. 

2. Let us demonstrate now that the structure (16.61) of spin operator is adequate to 
the considered model. Here we study a particular case that similar to one presented in 
Section 3, other cases can be solved in an analogous way. We take into consideration 
that the Hamiltonian H is given by structure the right-handed part of (16.21) and the 
spin operator we can find in the following form 



S 



( 







a 






S3 



o\ 

a 

Si J 



(6.7) 



We should check the main condition [S, H] =0. After substitution the specified 
expressions for H and S into this equation we obtain a system of linear equations 
for finding the coefficients sf. 



s 2 - si = 2p; 

s 4 - s 3 = 2p; 

S5 + s 6 = 0; 

m(si - s 3 ) + 2ps 5 = 0; 

m(s2 — S4) — 2psQ = 0. 



(6.8a) 
(6.86) 
(6.8c) 
(6.8d) 
(6.8e) 
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S 



\ 



(6.9) 



One of the solutions can be taken in the form s% = s 3 = —p, s 2 = s 4 = p, s 5 = s 6 = 0, 
and the corresponding spin operator can be presented as follows 

(-p a 0\ 

a+ p o o 

— p a 
a + pj 

We see that obtained spin operator is similar to one used in Section 3. 

Remark 1. The symmetry of the Hamiltonian can be violated as described in 
Sections 1-4. There we deal with the Hamiltonian that differs from (16.21 right) by the 
term -B7 5 , and B is a constant. In that case, the spin operator can be constructed in 
the form 

0\ 




/ 



s 



a, 
* 










(6.10) 



and the spectrum depends on spin quantum number similar to eq. (15. lip . Note that 
the spin operator (13. ip is similar to that in eq. (I6.10p . 

Remark 2. The theorem can be generalized to the Hamiltonians of more complicated 
structures namely 



H 



All 





h 3 


a\ 





h 22 


a + 




his 


a 


h3 







h 2 4 








or 



H 



All 


a 


h 3 


o\ 


a+ 


h 2 2 





h 2 4 







h 3 3 


a 


v° 


h 2i 


a+ 





(6.11) 







(6.12) 



The equation for the spectrum for this model (e.g., for the left side of eq. (16. 2p ) is 

h n -E h 13 /_(H) 

h 22 -E / + ((n)) h 24 

hi 3 /-((n)) h 33 -E 

U((n)) h 2A h M -E 

However, to solve this equation is not a simple task. 

Remark 3. It is not obvious from the matrix form of the Hamiltonian (15. 8 p , that it is 
a self-conjugate operator. To prove this, the properties of "increasing" and "decreasing" 
operators should be considered more carefully. A sequenced collection of functions 
Fj = y/\e\B£ l s (y^-r 2 ^ constitute a basis in the Hilbert space with scalar product 
defined as 

/•oo 

(F l s „F l s „)= / F l s ,F l s „rdr. (6.13) 
Jo 

So that, we get for each s and / 

(F l -\R-F l s ) = -(F l s ,R+F 1 - 1 ). (6.14) 
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Hence for operators ( 15. 7p we obtain {R~)* = —R + and (R + )* = —R~, where symbol * 
implies Hermitian conjugation of operators. From the matrix form of the Hamiltonian 
( 15. 8 j) we see that H* = H. So that we can use the theorem for this case. 

Note that the increasing and decreasing operators can be used effectively in this 
problem because we need in only one basis C l s for constructing the spinor ( I5.10p . In more 
complicated models, we seem should take two or even more different basis functions for 
constructing a solution. 



7. Spin coefficients 



We go back now to the Hamiltonian (I3.2p . In order to find its eigenvectors, we consider 
the following chain of transformations 

(-mT° + Gn - p )d + mC 3 = 0, 



(7.1) 



/ rro Gn Gn 
(mT - — +Po-—)C 1 



Po 




~ _ Gn 



mC 3 , 



Gn 



(mT° - 

Gn\2 
2 / 



(Po 



exH + — t^-Ci 



~ _ Gn 
Po ~ 



~ _ Gn 

Po ~ 



1- 



Hence, we obtain 



Ci 



mT° - % 



1- 



~ _ Gn 

Po ~ 



-A, C 3 = eJl + 



~ _ Gn 

Po ~ 



-A. 



In the same way we obtain the following expressions for Ci and C4: 

C 2 = 



~ Gn 

Po - — 



-B, C 4 = eWl + 



Gn 



~ Gn 

Po - — 



B. 



C 8 , 
(7.2) 

(7.3) 

(7.4) 



where A and B are new constants. 

Getting of the relation between G\ and C2 (and between C 3 and C4) is similar. 
Namely, we should use the equations 

(p 3 - mT°)C li3 + v / 2e 0J B(/ + s)C 2A = 0. (7.5) 

We finally get 



.4 



1 + -?LC, B 



1 - 



P3 



(7.6) 



with the only coefficient C which has to be defined from the normalization condition 
C\ + Cf + Cf + Cf = 1. We obtain G —\. Finally, we obtain the wave function: 

[C 1 C 1 - 1 ( e -fr 2 ) e^-D^ 



*(t,a;,y,«) = e" 



-ip t 



1 



2tt 



iC 2 4 (^fr 2 ) e 1 
C34- 1 (f f 2) e i(i-i)^ 

V iC 4 4(^r 2 )e^ y 



(7.7) 
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where 

1 / mT° - £p / p~7~ s ' I mT° - ^ / ^~ 

Cl = 2V 1_ po-^ V 1 + ^' C2= 2V 1_ po-^ V 1 "^' (7 ' 8) 

and L is a normalizing factor. 

Equations ()7.7p -( 17^|) represent the exact solution of (I2.4q|) with the Hamiltonian 
(12. 5p that describes the electron moving in matter and magnetic field. Note that in the 
case n — 0, these formulas are reduced to well-known solutions for the electron wave 
functions in a constant homogeneous magnetic field [Tj. 



8. Application 

It is important to point out, that the obtained exact solution for the electron motion 
in matter and magnetic field can be used as the first approximation of description of 
particles moving in an external environment of more complicated configuration. As an 
example, we demonstrate how the problem of an electron (or another charged particle) 
motion in a rotating matter with magnetic field can be solved. This problem is of 
interest in different astrophysical contexts. 

If the angular velocity u is small compare to the magnetic field, we can calculate 
the spectrum using a standard perturbation theory with the small parameter <C 1. 
If, for example, we choose for the matter density, angular velocity and magnetic field 
the values peculiar for a rotating neutron star (n = 10 37 sm~ 3 = 7.72 • 10 22 (eV) 3 , u = 
2tt ■ lO 3 *- 1 = 2tt • 0.66 • l(T 12 eV, B = 10 w Gs = 7 ■ 10 8 (eV) 2 ), then the parameter is 
really small, 

Gnu o n 

6.3-10~ 20 <l. 8.1 



e B 

Now we can take the spectrum and wave functions found above (i.e., without rotation) 
as the lowest order of perturbation series and find the correction term. 

We consider the particular case of constant magnetic field and rotating uniform 
matter so that the electromagnetic field and effective matter potential are given by 

^ = (0,-^,^,0), r = -Gn(l,-uy,ux,0), (8.2) 

where G = Let us rewrite eq. ( 12. 2 p in the Hamiltonian form 

i^*{x) = H*{x), (8.3) 

Gn 

H = 7° 7 (p + e A) + m 7 ° + — (1 - 4 sin 2 9 W + 7 5 ) + (8.4) 
+ ^T7V(1 - 4 sin 2 9 W + l b )ujy - ^7°7 2 (1 - 4 sin 2 9 W + -f 5 )ux, 



Relativistic particles in magnetic field and matter 12 

with A = (— 0). Using again the chiral representation of the 7-matrices we obtain 
the Hamiltonian in block-matrix form: 



H 



—a(p + e A) + Gn m 

m cr(p + e A) 

-2Gnsm 2 6 w -Gnu(^ ° lV ~ ° 2X °V (8.5) 

It is obvious from (18. 5p . that H = H + Hi, where Hi is the last term of (18. 5 j) and 
takes the following form in the polar coordinates 



Hi 



( 


— ipre 1(1 


' 


o\ 


ipre 1 ^ 























I 








V 



p = Gnu. (8.6) 



So, we get for the first correction to the energy spectrum (14.21) of the electron 

00 

Ap% = J ** N Hi* N dV = ^2e~BCiC 2 2 -^- J C^O^n-M) (8-7) 



where Ci and C2 are the coefficients from eq. (17.81) . The integral can be calculated 
explicitly, and finally we get 

2Gnu 



Ap$ = y f 2e^BC l C 2 — V^V 

e -D 



= 2GnujC x C 2 J^- = 2GnC l C 2 uj(R). (8.8) 
V e B 

This shift of levels in the energy spectrum depending on the energy quantum 
number iV = 0,1,2... leads to a corresponding shift in a frequency of synchrotron 
radiation of electron inside of dense rotating matter, that can be registered. 



9. Conclusion 



In this paper we found a class of exact solutions of the modified Dirac equation 
which describes the charged leptons propagating in uniform matter and strong constant 
magnetic field. We also pointed out, how this approach can be generalized to a given 
class of Dirac Hamiltonians. Obtained solution for the particular case of the electron 
motion in a rotating neutron star with account for matter and magnetic field effects 
can be used as the first approximation in more complicated models. All of these 
considerations can be useful for astrophysical applications. 



References 



[1] Sokolov A A and Ternov I M 1968 Synchrotron radiation (Oxford: Pergamon Press) 



Relativistic particles in magnetic field and matter 



13 



[2] Ritus V 1978 Issues in intense-field quantum electrodinamics ed by V Ginzburg (Moscow: Proc. 

of Lebedev Phys.Inst. Ill 5) 
Nikishov A 1978 Issues in intense-field quantum electrodinamics ed by V Ginzburg (Moscow: Proc. 

of Lebedev Phys.Inst. Ill 152) 
[3] Furry W 1951 Phys.Rev. 81 115 

[4] Studenikin A 2006 Ann.Fond. Louis de Broglie 31 no. 2-3 289 

Studenikin A 2008 J.Phys.A: Math.Theor. 41 164047 

Studenikin A 2006 J.Phys.A: Math.Theor. 39 6769 

[5] Studenikin A and Ternov A 2005 Phys. Lett. B 608 107 [Preprint [hep^ph/0410297|) 

[6] Grigoriev A, Studenikin A and Ternov A 2005 Phys. Lett. B 622 199 (Prepnnf |hep-ph/0502231 1 

Grigoriev A, Studenikin A and Ternov A 2005 Grav. & Cosm. 11 132 

Grigoriev A, Studenikin A and Ternov A 2006 Phys. Atom. Nucl. 69 1940 
[7] Studenikin A 2006 Annates de la Fondation Louis de Broglie 31 289 

Grigoriev A, Studenikin A, Ternov A, Trofimov I, Shinkevitch S 2007 
Russian Physics Journal 50 596 

Grigoriev A, Shinkevitch S, Studenikin A, Ternov A, Trofimov I 2008 Grav. & Cosm. 14 248 
[8] Grigoriev A, Savochkin A, Studenikin A 2007 Russian Physics Journal 50 845 
[9] Balantsev I, Popov Yu, Studenikin A 2009 II Nuovo Cimento 32 C 53 



